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Abstract 



It is highly plausible that the region of space-time far from an isolated gravitating body is, 
in some sense, asympt otically Minkowsk i an. However theo retical studies of t he fu ll nonlin- 
ear theory, initiated by [Bondi et all (119621) . ISachsl (119621) and iNewman & Until (119621) . rely on 
careful, clever, a-priori choices of chart (and tetrad) and so are not readily accessible to the 
numerical relativist, who chooses her/his chart on the basis of quite different grounds. 

This paper seeks to close this gap. Starting from data available in a typical numerical 
evolution, we construct a chart and tetrad which is, asymptotically, sufficiently close to the 
theoretical ones, so that the key concepts of Bondi news function, Bondi mass and its rate 
of decrease can be estimated. In particular these esimates can be expressed in the numerical 
relativist's chart as numerical relativity recipes. 
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1 Introduction and motivation 



The two threads which underpin this study of asymptotic flatness are theoretical and numerical rel- 
ativity. We start by reviewing the former. It is widely believed that the region of space-time far from 
an isolated gravitatin g body is, in some s ense, asympt otically Minkowskian. Already in 1962 Bondi 
and his coworkers (|Bondi et all 1 19621) . (ISachslll962l) developed asymptotic expansions for the so- 
lution of the full nonlinear vacuum field equations, leading to a r igorous concept of gravitational 
radiation in the far field. Soon after wards Newman & Until (119621) produced an alternative version 
using the NP null tetrad formalism (fNewman & Penrose! . 1 1962|) . A key ingredient in this and later 
work was the careful choice of a suitable coordinate chart involving a "retarded time" coordinate 
u. Both groups introduced a (u, r, 9, <p) chart where 9 and <p were spherical polar coordinates. They 
both developed as ymptotic expan sions as r — ► oo holding the other coordinates fixed. Subsequent 
work by Penrose ([Penrose! . Il963|) showed that by a process of "conformal compactification", in- 
finity could be adjoined to the space-time manifold and then treated by standard methods. Most 
modern theoretical treatments use the Penrose conformal approach, and the "old-fashioned" chart- 
based approach has fallen out of fashion. However it is closer to what most numerical relativists 
are calculating, and for this reason we shall use it here. 

The two groups used different charts and dependent variables. iBondi et al.1 (|1962|) chose a 
(u, r, 9, 0) chart where r was an area coordinate and 9 and <p were standard spheri cal polar coordi- 
nates (see below). Their primary dependent variables were the metric components. [Newman & Unti 
dl962h produced an alternative version using a null "retarded time" coordinate u. Then the null 
vector l a = g ab u^ is geodesic and their coordinate r was chosen to be an affine parameter for the 
integral curves of l a along which the other three coordinates were fixed. Their primary dependent 
variables were the tetrad connection components and the tetrad components of the Weyl curvature 
tensor. These ten independent Weyl tensor components are usually described by five complex scalar 
funct ions where n = 0, 1, . . . , 4. (For a covariant physical interpretation of the Weyl tensor see 
e.g.. lSzekeresldl965lU 

Both groups wer e considering th e limit r — > oo with the other coordinates fixed. In the Penrose 
geometrical picture ([Penrose! . 1 19631) this region is called future null infinity. Of course both groups 
could have considered an "advanced time" coordinate v where the corresponding limit is past null 
infinity. The \l/ n in that case have similar properties and interpretation to ^i- n near future null 
infinity. 

If one wants to consider an isolated system with no extrinsic incoming radiation, then, as ex- 
plained below, the natural place to impose this is past null infinity. How ever both groups looked 
for a condition to be imposed near future null infinity. IBondi et all (|1962|) introduced an "outgo- 
ing radiation condition" which required the vanishing of certain terms in the asymptotic expan- 
sion of two of the metr ic components. This condition will be stated more precisely in section [3] 
Newman & Unti dl962l) made a "peeling assumption": near future null infinity \& = 0(r 5 ), and 
with this assumption they were able to demonstrate a so-called "peeling theorem": = Q(r n ~ 5 ) . 
Then ^ 4 = 0(r _1 ) is interpreted as the leading term in the outgoing radiation. In the (IBondi et al. . 



1962) picture the equivalent role is taken by the "Bondi news function" built from first derivatives 
of metric components. For a comparison of the conditions in the tw o schemes, showing th at the 
outgoing radiation condition implies the peeling assumption see e.g., IValiente-Kroonl (|1999|) . and 
section [51 



2 



By reversing the direction of time, swapping advanced time for retarded time, one could carry 
out an almost identical study near past null infinity. There, assuming the analogous peeling condi- 
tion, \T/ = 0(r _1 ) is to be interpreted as the leading term in the extrinsic incoming radiation, and 
a natural "no incoming radiation condition" near past null infinity would be ^4 = 0(r~ 5 ). 

It is important to realise that the "outgoing radiation condition" or "peeling assumption" does 
not preclude the presence of incoming radiation near future null infinity. Even within linearized 
theory the "peeling theorem" allows modest amounts of incoming radiation (Deadman & Stewart, 
paper in preparation). 

We turn now to numerical relativity where researchers have expended considerable effort on the 
(numerical) evolution of asymptotically flat space-times. A minority of researchers have adopted 
the Penrose conformal approach, but most have chosen to evolve the space-time as far out (both 
in space and time) as is feasible, using the traditional approach. Then some matching process 
is required to interpret their numerical data in the Bondi or NP pictures. This, the goal of this 
paper, turns out to be far from trivial. The choice of a coordinate chart is an intrinsic part of 
the numerical evolution and the final data is available only in this chosen chart. Each numerical 
relativity group has its own favoured chart or charts and they usually bear little resemblance to 
the Bondi or NP ones. Furthermore this data does not contain complete information because the 
inevitable occurrence of numerical errors will corrupt the values of higher derivatives — from it one 
can construct reliably only a few leading terms in asymptotic expansions^. 

The usual approach adopted by numerical relativists is to argue that far from the isolated source 
the gravitational field i s weak, and s o line arized theory can be used to match the numerical and the 



Bondi or NU pictures. iBondi et al.l (11962 ) argued strongly against such an approximation pointing 
out the fundamental nonlinearity of general relativity. Even if plausible arguments in its favour 
could be found, linearization carries its own difficulties. The first is that in a non-compactified 
space-time the matching process is a global one. Further given a space-time, the choice of a sim- 
pler second space-time of which the first can be considered a linearized perturbation, is not unam- 
biguous. Even if such a choice could be justified the transformation between the charts in the two 
space-times would not, in general, be smooth. 

As a concrete example illustrating these points consider the well-known Schwarzschild metric 
in the standard (t, r, 9, 4>) chart 

g s ab = diag(F, -F-\ -r 2 , -r 2 sin 2 9), (1) 

where F = 1 — 2M/r. In the region where r ^> M this might appear to be a small perturbation of 
Minkowski space-time with metric 

gf b = diag(l, -1, -r 2 , -r 2 sin 2 9), (2) 

but this is deceptive. Consider the scalar wave equation g ab ^ ]a b = on the two space-times. We 
would measure outgoing radiation at future null infinity by taking the limit r — > 00 holding u 



'Consider an asymptotic expansion 

f(r) =fo + fir' 1 + f 2 r- 2 +... 

as r — > 00. If we interpret this as the first few terms in the Taylor series for /(g) about q = where q = r _1 , then the 
/„ are, up to numerical factors, the (/-derivatives of / evaluated at q = 0. 
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constant, where u is a retarded time coordinate. Two standard choices for u are 



u M = t - r, u s = t - r*, (3) 

where 



F~ L dr = r + 2M log \r/2M - 1| + const. (4) 
Thus 

u M = u s + 2M\og\r/2M - 1\+ const. (5) 

The Schwarzschild null infinity is given by r — > oo holding constant, which implies u M — ► oo, 
known as future timelike infinity for the Minkowski space-time. Equivalently the Minkowski null 
infinity involves taking the limit with u M constant which corresponds to r — > oo with u s — > — oo, 
known as spacelike infinity for the Schwarzschild spacetime. Thus the limits in the two charts are 
different. This happens because of the global nature of the limiting process. 

In order to achieve comparable limiting processes we need to redefine the two charts. Here 
both space-times are static and so it is simplest to retain the t-coordinate. Suppose we invert (for 



> M) the relation [(4) , r* = r*(r) giving r = r(r*) and introduce a new chart (t, r*, 9, <ft). Then 



the Schwarzschild line element [(T)1 becomes 

g s ab = diag(F, —F, -r\ -r 2 sin 2 9). (6) 
Using the same chart the Minkowski line element is 

g% = diag(l, -1, -r 2 , -r 2 sin 2 9). (7) 



Now the two metrics |(6)| and [(7)] are not only small perturbations of each other (for large r*), but 
they share the same causal structure, u = t — r* in both cases. (There are of course many other 
ways of doing this, e.g., retain the r's and change the t's, which is the approach to be adopted in 
this paper.) Note also the appearance of logarithms, which means that the transformations are not 
smooth. 

The purpose of this paper is to examine in more detail these issues from the point of view of 
the numerical relativist. In section [2] we state what information we believe is available in a typical 
numerical evolution, and we assume that this information is expressed in terms of a given chart 
X a = (T, R, 6,$) which is asymptotically Minkowskian. Section [3] addresses the construction 
of an approximate Bondi-like chart x a = (u, r, 9, <p) using this information. This circumvents the 
problem referred to above. We write down here the explicit fo rm of the Bondi et al. outgo ing radi- 



ation condition. We introduce a Newman- Penrose (NP) tetrad (|Newman & PenroseL 11962) adapted 



to the problem by [Newman & Unti dl962h in section |4l At leading order this is the usual NP tetrad 



for Minkowski space-time. At each order, r _1 , r~ 2 , . . . , there are 16 real coefficients describing 
the tetrad. However from section [3] we know that only 10 coefficients are needed to describe the 
metric. There are six coefficients which describe an infinitesimal Lorentz transformation at each 
order, and, for the moment, we do not make a particular choice for them. 

In section |5] we first obtain the asymptotic solution of the full nonlinear vacuum Einstein equa- 
tions. As we do so we fine-tune our chart and NP tetrad to make them closer to those of Bondi et 
al. and Newman-Unti. Once we have set the Ricci curvature, to the best of our abilities, to zero 
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we turn to the Weyl curvature described by the Weyl scalars \I/ n referred to earlier. We find that 
= 0(r n ~ 5 ) for n = 4, 3, 2, 1, but \l>o = 0(r~ 4 ), which would appear to violate the NU peeling 
assumption. However using the information gleaned from solving the vacuum field equations and 
the fine-tuning of the chart and tetrad, we can show that the Bondi outgoing radiation condition 
implies the NU peeling condition so that the peeling theorem then holds. 

The bad news is that the leading terms in the Weyl scalars \& = 0(r~ 5 ) and = 0(r~ 4 ) 
cannot be estimated using the information we judge to be available from the information extracted 
in section [2] Although these scalars can be computed in linearized theory, that theory would appear 
to be an unreliable guide here near future null infinity. 

The good news is that we can compute the leading terms in $4 = 0(r _1 ), equivalent to the 
"Bondi news function" (and we can compute this scalar accurately within linearized theory). The 
same holds for ^3 = 0(r 2 ), which involves nonlinear terms, but these can be removed by the fine- 
tuning process. We can also compute ^2 = 0(r~ 3 ) which involves nonlinear terms in an essential 
way. This means that we can offer reliable estimates of the "Bondi mass" Mb(u) of the isolated 
systerrH and its rate of decrease dMs/du ^ 0, presumably due to the radiation of energy, both 
manifestly inaccessible to linearized theory. 

The final section [6] translates these results back into the X a chart of section [2] used by a typ- 
ical numerical relativist. From her/his standpoint there is no need to go through the elaborate 
construction of a theoretical chart and NP tetrad carried out in the intermediate sections. We offer 
"numerical relativity recipes" so that they can compute the key quantities referred to in the previous 
paragraph in their own preferred chart. 

The key ideas in this paper are at least forty years old, a nd one migh t ask why wer e these 



resul ts were not give n before? The nonlinear calculations of iBondi et al.1 (| 1962T) . ISachsl (I1962|) 



and Newman & Unti dl962h were made possible by careful, clever, a-priori choices of chart and 



tetrad. We have to start from more or less arbitrary choices and so the resulting expressions are 
horrendously complicated. In order to handle them accurately we have utilized a computer algebra 
system. We used Reduce, and our Reduce 3.8 scripts can be obtained by an email request to the 
authors. Our choice reflected our experience and knowledge of one particular computer algebra 
system, but we used no features not available in some other systems. 

We are grateful to Oliver Rinne for very useful discussions. One of us (ED) is grateful for 
financial support from the UK Engineering and Physical Sciences Research Council. 



2 The "Bondi mass" is of course the timelike component of a 4-vector and so frame dependent. But a numerical 
relativity evolution singles out a well-defined frame, and that is the one in which the mass is computed. 
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2 The numerical data 



Most numerical relativists would choose a quasi-spherical polar chart X a = (T, R, 6, $) for the 
numerical evolution of the space-time surrounding an isolated gravitational source. We could also 
define an associated quasi-cartesian chart Y a = (T, X, Y, Z) where 

X = .RsinOcos$, Y = Rsin 6 sin <3>, Z = RcosQ. 

We shall be interested in the limit R — > oo. As stated this limit is meaningless unless we 
specify the behaviour of the other three coordinates under the limiting process, and we shall rectify 
this omission shortly. It proves very convenient to introduce the notation 

O n = 0(R~ n ) as R -> oo. (8) 

Our fundamental assumption is that the space-time outside an isolated source is asymptotically 
Minkowskian, expressed by the idea that, as seen in the Y a chart, 

9* = Vab + gStRT 1 + gVRT 2 + 3 , (9) 

where r] ab = diag(l, —1,-1, —1) and the are supposed to remain constant during the limiting 
process. Transforming from the Minkowskian chart to the spherical polar one we find that the 
metric components in the X a chart look like 

#oo = 1 + hooR' 1 + kooR' 2 + O s , 

0oi = hoiR- 1 + k 01 R- 2 + 3 , 

go2 = h 2 + k 02 R~ 1 + 2 , 

903 = h 03 + k^R' 1 + 2 , 

9 u = -1 + h lx R~ l + k n R~ 2 + 3 , 

gu = h 12 + k 12 R^ + 2 , 

9i3 = hi 3 + k^R" 1 + 2 , 

922 = -R 2 + h 22 R + k 22 + Oi, 

923 = h^R + ha + Oi, 

933 = -R 2 sin 2 6 + h 33 R + k 33 + O x . 

Here the functions {h ab } and {k a b} are required to remain constant during the limiting process. 
We shall also need the asymptotic form of the inverse metric g ab which is readily obtained from 
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the relation g ac g c b = <^V We find 



g 00 = 1 + h 00 R -l + k 00 R -2 + 

g <* = h m R~ l + k 01 R~ 2 + 3 , 

g™ = h 02 R -2 + ^02^-3 + 

g 11 = -1 + h 11 R~ 1 + k n R~ 2 + 3 , 

= /, 12 iT 2 + k 12 R~ 3 + 4 , 
3 13 = h 13 R~ 2 + k 13 R- 3 + 4 , 

^22 = _^-2 + ^22^-3 + £.22^-4 + ^ 

#* = h 23 R~ 3 + k 23 R' 4 + 5 , 

g^ = -R- 2 CSC 2 6 + h 33 R~ 3 + k 33 R~ 4 + 5 . 



Explicit formulae for h ab and k ab are given by equations |( A 1 )| and |( A2)| in appendix lAl At this level 
of approximation 

g ac g cb = 5\ + o 3 . 

A numerical evolution in which the dependent variables include both g ab and g abc , usually 
called a "first order formulation", should produce accurate values for h ab and its first derivatives, 
and for k ab . Otherwise we assume that these variables are available for discrete ranges of T, and 
$ so that the corresponding derivatives can be estimated. 
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3 The Bondi chart 



Most of the theoretical work which has been done on outgoing gravitational radiation involves a 



"Bondi chart" (u,r,9,4>) i n which u is a retarded time coordinate, see e.g., iBondi et al.l (|1962f) . 



Newman & Penrose! d!962h . iNewman & Until (119621) 



Here we take the viewpoint that the X a = (T, R, 6, $) chart introduced in section [2] is the 
fundamental one in which, ultimately, all numerical calculations will be performed. Starting from 
this chart we need to construct an x a = (it, r, 9, <p) one which has all the essential features of a 
Bondi chart and we start by studying the function it(T, R, 0, $). 

Because it is a null coordinate it has to satisfy the relativistic eikonal equation 



g a \ a u, b = 0. (12) 

This is a well-known nonlinear equation with four independent variables which is exceedingly 
difficult to solve with any generality. (Even the restriction of |(12)| to Minkowski space-time leads 
to the surprisingly rich structure of light ray caustics.) Note that there is a "gauge freedom" — if u 
is a solution then so is U (u) for any differentiable function U. 

The standard approach is to specify u on a spacelike hypersurface in space-time, and then 
existence and local uniqueness of u is guaranteed by standard theorems. The standard approach 
is of little utility in this context, for no obvious choice of data suggests itself, and so we adopt a 
different approach. 

Consider first the special case of a Minkowski space-time, where |(T2)1 can be rewritten as 

(u T ) 2 - (u, R ) 2 = R- 2 [(u e ) 2 + esc 2 6(ii, $ ) 2 ] = 2 . (13) 

Suppose we look for spherically symmetric solutions u = u(T, R). Setting uj = u.r/u^t in |(13)| 
we find uj 2 = 1. Using the gauge freedom mentioned earlier we may impose it T = 1 to find 

du = dT ± dR, 

which implies 

u = T ±R + const. 
T — R is called retarded time and T + R is called advanced time. 

Although the special case appears trivial it is the key to the general one. Within this section 



only let indices i, j range over 0, 1 and let indices /, J range over 2, 3. Perusal of the display [(TT) 
shows that g % i is Oq while both g' lJ and g IJ are Oi- Thus the eikonal equation takes the form 

g^UiUj = 2 , (14) 

which should be compared with |(13)| above. As boundary conditions (as R — > oo) we impose 

u, T = l + O u v = °i- O 5 ) 
This means that the eikonal equation takes the form 

g ij UiUj = 3 , (16) 
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which we can write as a quadratic equation for uj = U r/u t> an d choosing the sign appropriate for 
retarded time we find the solution 



where 
and 



u R = -(1 + 2 mi /R + 2m 2 /R 2 )u T + 3 , (17) 
m 1 = -\{h 00 + 2h 01 + h n ), (18) 



m 2 = - yg [4/coo + 8fc i + 4fcn + (h 00 - h u ) 2 - 

4(/i o + V) 2 + 4(/i 02 + h 12 f + 4(/i 03 + h 13 ) 2 esc 2 



(19) 



Thus 



du = u jT dT - u T (l + 2mi/i2 + 2m 2 /R 2 ) dR + 3 . (20) 
We leave some freedom in it by setting 

^ = 1 + l + f + - (21) 

where q± and ^2 are i?- independent functions. At the moment they are arbitrary. The require- 
ment that the vacuum Einstein equations hold then determines inter alia q\, see section 6. In our 
calculation q 2 is not used directly. 

We need next to specify a radial coordinate r = r(T, R, 0, $), and the simplest choice is r = R. 
This has the great practical advantage that O n = 0(R~ n ) = 0(r~ n ). It could be argued that our 
choice of r is n either the Bondi area co ordinate nor an affine parameter along the outgoing null rays 
as favoured by iNewman & Until (| 19621) . However since both of those approaches are known to be 
essentially equivalent, it would seem that the discussion is not sensitive to the precise choice of r. 

Then |(20)| implies 

dT = ((1 - q x /r - (q 2 - qi 2 )/r 2 ) du + (1 + 2m x /r + 2m 2 /r 2 ) dr + 3 , (22) 

and so 

rap, . 1 _*_£ S £ 1 ! +03 , rm =0i (23) 

\ou J r r z \ou J r 

(&T\ 2m 1 2m 2 „ fdR\ 

TP =l + ^ + ^ + 3 , U- =1. (24) 

Finally we consider the choice of angular coordinates 9 = 9(T, R, 6, $) and = 0(T, R, 9, $). 
We shall require 9 = + 0\ and = $ + Oi, and so the relations, being close to the identity, are 
invertible. It is more convenient to posit 

Q = 9 + ^ + ^ + 3 , $ = 0+^ + ^| + O 3 , (25) 

where the functions yj and zj do not depend on r but are otherwise arbitrary. Equations |(25)| are 
certainly consistent with the boundary conditions |(15)| 
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We can now specify the limiting process as r — > oo holding u, 9 and </> constant. Thus we are 
regarding m n , q n , yj, Zj, {h ab } and {k ab } as functions of u, 9 and (p. 

Because we know the Jacobian (dX a /dx b ) we can write down the metric components in the 



(u, r, 9, 0) chart 



5-00 = 1 + a 00 /r + b 00 /r 2 + O 



3) 



0qi = 1 + a i/r + boi/r + 3 , 

902 = -ry2 lU + a 02 + W r + °2, 

903 = -rz 2iU sin 2 9 + a 03 + b 03 /r + 2 , 
gn = a n /r + b n /r 2 + O s , 

9X2 = an + b u /r + 2 , 
9i3 = Oi3 + b 13 /r + 2 , 

922 = ~r 2 + a 22 r + b 22 + O h 

923 = a 23 r + b 23 + O u 

933 = -r 2 sin 2 9 + a 33 r + b 33 + O x . 

Two points should be noted here. Firstly the leading terms in g 02 and g 03 , if non-zero, would violate 
our notion of an asymptotically Minkowskian space-time, for they are not present in the standard 
Minkowski line element. Thus we need to impose the conditions or "constraints" 

V2.u = Z2,u = 0. (27) 

Explicit formulae for a rnn in terms of the h mn , g ls m 1 and y 2 , z 2 (after imposing |(27)[ ) are given as 
|(A3)| in appendix lAl We could also give explicit formulae for the b mn in terms of the h mn , k mn , q n , 
m n , y n and 2;„but they are rather lengthy, and are most easily generated using a computer algebra 
package. 

Next recall that the u-coordinate was constructed as a solution of the eikonal equation |(16)[ 
Thus, as seen in the (u, r, 9, cp) chart g 00 = 3 which implies g u = 0\ and so a n = 0. One may 
verify this directly by comparing the explicit expression for a n given in |(A3)| with |( 1 8)[ We will 
show later that by making a suitable choice for y 2 and y 3 , we can achieve bn = so that gn = 3 
as expected. 

We now have s ufficient notation available to write down the "outgoing radiation condition" of 



Bondiet all (119621) as 

a 33 = ~a 22 sin 2 9, b 33 = b 22 sin 2 9, b 23 = 0, (28) 
which we shall invoke later. 
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4 The NP tetrad 



Since most recent studies of gravitational radiation use a NP null tetrad (Newman & P enrose!. 1 1962|) 
we need to introduce one. The basics of tetrad formalisms are du e tolSchoutenl (H 954). Many text- 
books contain more readable, but often succinct accounts, and Chandrasekh ail dl983b chapter 1, 
section 7, is a good pedagogic compromise. With small, but necessary, changes in notation this 
is summarised in appendix El The special i satio n of this approach to the original NP formalism 
has been given by ICampbell & Wainwrightl (|1977|) . It turns out that the calculations that we need 
to do with it become surprisingly intricate, and so are most conveniently handled using a com- 
puter algebra system. One of th e authors implemented the appendix |B] formalism, the other the 
(|Campbell & Wainwrightl. 1 1977|) one. Both of course gave the same results, a useful guard against 
programming errors. Results here are reported for the former. 

We use a tetrad of vectors e a a and the dual tetrad of co vectors e a a . (The tetrad indices are greek 
characters and always occur first.) Tetrad indices are lowered and raised using e a p and e a/3 where 





/0 1 





\ 


c — — 

e«/3 — e — 


1 









-1 




\ 


-1 


o / 


In NP notation we have 








e° = I e 1 = n 


e\ = 


m a , 




e a = n a , ei a = Z a , 


e 2 a = 


—m a 


, e 3 ' 



m„ 



-m 



(29) 



We shall require that, to leading order, l a = w , 
Setting s = 2~ 1 / 2 we write the tetrads a^U 



e° a = (1 + coo/ r + doo/r 2 + 3 , c 01 /r + d 0l /r 2 + 3 , 

c 2 + d 02 /r + 2 , c 03 + d 03 /r + 2 ), 
e l a = {\ + cio/r + d w /r 2 + 3 , 1 + c n /r + d n /r 2 + 3 , 

C12 + du/r + 2 , ci 3 + d 13 /r + 2 ), 
e 2 a = (c 20 /r + d 20 /r 2 + 3 , c 21 /r + d 21 /r 2 + 3 , 

- sr + c 22 + d 02 /r + 2 , isr sin 9 + c 23 + d 23 /r + 2 ), 
e 3 a = {c 30 /r + d 30 /r 2 + 3 , 1 + c 31 /r + d 31 /r 2 + 3 , 

- sr + c 32 + d 32 /r + 2 , -isr sin 9 + c 33 + d 33 /r + 2 ), 

3 In our calculations we actually included one extra term in each of the asymptotic expansions below. E.g., the first 
component of e° a was written as 

e% = 1 + Wr + d 0Q /r 2 + j 0Q /r 3 + 4 . 

These "junk" terms show up in our expressions for the connection and curvature components. In any expression where 
a junk term occurs we regard all terms of that (and any higher order) as being junk, and not computable form the data 
described in section|2] 
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and 



e a Q = (1 + c 00 /r + d 00 /r 2 + 3 , -\ + c 01 /r + d 01 /r 2 + 3 , 

c 02 /r 2 + d 02 /r 3 + c 03 /r 2 + rf 03 /r 3 + 

el = (c 10 /r + d 10 /r 2 + 3 , 1 + c u /r + d u /r 2 + 3 , 
c 12 /r 2 + d 12 /r 3 + 4 , c 13 /r 2 + d 13 /r 3 + O4), 

e « = ( C 20 /r + d 20 /r 2 + c 21 /r + rf 21 /r 2 + 03; (3D 

- s/r + c 22 /r 2 + d 02 /r 3 + 4 , -is esc 6/r + c 23 /r 2 + d 23 /r 3 + 4 ), 
= (c 30 /r + d 30 /r 2 + 3 , 1 + c 31 /r + d 31 /r 2 + 3 , 

- s/r + c 32 /r 2 + d 32 /r 3 + 4 , is esc 5/r + c 33 /r 2 + d 33 /r 3 + O4). 

Each tetrad contains, at each order, 32 real coefficients. For although c 2n and C3 n are complex, 
= c 3n etc. The relation e^ a e v a = 5^ u allows one to determine the c mn in terms of the c mn and 
the d mn in terms of c mn and d mn , reducing the number of unknowns, at each order, from 32 to 16. 
The first set of these is given as equation |(A4)| in appendix lAl. The second set is rather lengthy and 
best generated using a computer algebra package. 

We also have the relation e^ u e^ a e v b = 9ab, and this enables us to determine a mn in terms of c mn . 
Note that there are 10 a mn and 16 real c mn . Given the tetrad, the metric is uniquely determined. 
But for a given metric there is a 6-parameter set of tetrads which give rise to it. They are of 
course Lorentz transformations of each other and the Lorentz group has 6 arbitrary parameters. We 
introduce 6 arbitrary first order Lorentz parameters a m (u, 9, (j>) and can determine the c mn in terms 
of the a mn and a m . There are many different ways of doing this, and one is written down explicitly 
as |(A5)| in appendix lAl We can of course write down the d mn in terms of the a mn , b mn , a m and 
extra second order Lorentz parameters (3 m , but they are rather lengthy and are best generated by 
computer algebra. 
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5 The Curvature Tensors 



We now use the tetrads developed in section |4] to evaluate the Ricci and Weyl curvature tensors 
using the algorithm outlined in appendix [Bj At each stage we convert all instances of c mn and 
d mn to instances of c mn and d mn using |(A4)| and its second order analogue. Then we convert all 
instances of c mn and d mn to instances of the metric coefficients a mn and b mn using |(A5)| and its 
second order analogue. These conversions are implicit and will not be mentioned explicitly again. 

We start by looking at the Ricci tensor component R n = R ab ei a ei b = R ab l a l b . We find 

= R U = a lhu /r 2 + 3 . (32) 

We chose our chart to ensure that u was approximately a null coordinate or equivalently g 00 = 3 
or equivalently g u = 0\. This means that we have to enforce a n = 0, and so the leading term in 
Rn vanishes. We look next at 

= Rot = RabeoW = R ab n a l b = -|a U|UU /r + 2 . (33) 

Again the leading order term vanishes automatically. Next consider 

= R 12 + R 13 = -2sa 12iM /r 2 + 3 , = R 12 - R 13 = -2zsa 13 , u esc 9/r 2 + 3 , (34) 
where s = 2~ 1 ^ 2 . We deduce that 

«12,« = an,u — 0. (35) 

Further we can compute 

= R02 + R03 = sa 12tUU /r + 2 , = R 02 - Ro 3 = isa 13iU esc 9/r + 2 , (36) 
If we inspect R 23 and use |(35)| we find 

= R 23 = \{a 22>u + a 33:U esc 2 6)/r 2 + 3} (37) 

and we deduce that 

022,« + a 3 3,« esc 2 6 = 0. (38) 

Further we may compute 

#oo = -|(a22,n« + a 33tUU esc 2 9)/r + 2 , (39) 

and we see immediately from [(38)] that the leading term vanishes, and so furnishes no new informa- 
tion. Finally inspection of the leading 2 terms in R 22 ± R 33 reveals that they vanish automatically 
because of an = 0, |(35)| and |(38)[ Thus we have 

R 22 + R 23 = 3 , R 22 — R 23 = 3 . (40) 

We have found, so far, that the conditions an = 0, |(35)| and |(38)| imply that R 00 , R i, R 02 and R 03 
are 2 while the other components are 3 . 

At this point we need to examine [(3 5 ) | more closely. Using [(A3)] we have 

(h 2 + h 12 + y 2 ) >u = (h 03 + h 13 + z 2 sin 2 0) iU = 0. (41) 
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Now we know that the functions y 2 and z 2 are arbitrary, apart form the constraints \(27)\ and so |(41)| 
implies 

(h 02 + h 12 ), u = (h 03 + h 13 ) >u = 0. (42) 
We may therefore choose, consistent with the constraints |(27)[ 

V2 = -O02 + h 12 ), z 2 = -(h 03 + h 13 ) esc 2 9, (43) 

which, using |(A3)1 sets 

ai 2 = ai 3 = 0. (44) 

The choice [(43)1 has an added advantage that if we now express bn in terms of the h a b, k a b, mi 
and m 2 , y 2 and z 2 we find that bn = 0, so that gn = 3 . At the same time we can examine b\ 2 and 
6 13 which are linear in y 3 and z 3 respectively. By choosing y 3 and z 3 appropriately we may arrange 

&12 = b n = 0. 

This is a convenient point at which to examine the choice of a specific Lo rentz transformation, 
in our tetrad, exemplified at leading order by the parameters a n , see e.g., |(A5)] iNewman & Penrose 



(|1962l) chose u to be a null coordinate, g ab u^ a Ub = and the covector l a = u >a . For a symmetric 
connection it follows easily that l a ] t > l b = 0. Even if we set l a = f(x c )u A we find that l a $l b is 
proportional to l a so that we still have a null geodesic, albeit not necessarily an affinely parametrised 
one. The rest of this paragrap h relie s on some details of NP formalism which can be checked swiftly 



using appendix B of lstewartl (1993). Within NP formalism 

l a . b l b = ( e + z)l a - Rm a - nrrf ; (45) 

where k and e are NP spin coefficients defined below. Now our coordinate u is only approximately 
null, and our covector l a is only approximately its gradient. Here k = m a l b l a $ = 7131 (using the 
notation of appendix |B]) turns out to be O3. However if we choose 04 = a$ = (3^ = (3$ = 0, and 
anticipate a 01 = (see next paragraph), we obtain n = 4 . Also e + e = n a l b l a . b = 7 n is 2 but 
if we choose «i = we find e + e is 3 . We also found r = 7i 30 = (a 2 + ia 3 )/r 2 + 3 or r = 3 
if we impose a 2 = a 3 = 0. At this stage we also choose a 6 = (3 e = for reasons given below. 

Now we need to examine each of the remainder (next order) terms in |(32)[ |(33)[ |(34)[ |(36)[ 
|(37)[|(39)1 and |(40)j For example we now find that the 3 terms in R n vanish if and only if we set 
a i = 0, and then Rn = O4. This also implies that the 2 terms in i? i vanish, so that i? i = 3 . 
We already established that Ri 2 ± i?i 3 = O3. Setting the leading order terms to zero furnishes 
expressions for ao 2 and a 03 which we use for subsequent simplifications. Now R\ 2 ± Rn = O4. 
We find then that our previous estimate [(36)| refines to -R02 ± -R03 = O3. We also need to refine our 
estimate [(37)| to R 23 = O4. We find R 22 ± R 23 = O3, where both O3 terms deliver the same relation 
relating the n-derivatives of a 22 , a 23 , a 33 and 622, b 23 , b 33 which we save for later use. In deriving 
this result we had to choose = (3§ = and to impose the Bondi outgoing radiation condition 
|(28)[ Then R 22 ± R 23 = O4. Next we reexamine |(39)[ The leading 2 term gives us an expression 
for a 00tU which we store for later use. Now R 00 = 3 . We have found, so far, that R 00 , R i, R 02 
and R 03 are 2 while the other components are 3 . 

We now try to repeat the procedure of the previous paragraph. However we find that the O4 
contribution in R u contains "junk" terms, i.e., terms which involve the third order metric com- 
ponents which we have not been including; see the footnote in section HI Thus we can obtain no 
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further information from the vacuum field equation R n = 0. Similarly we find that the 3 terms in 
Roi contain junk as do the O4 terms in R\ 2 ± -R13. The same applies to the 3 terms in R 02 ± R03, 
the O4 terms in R 23 and R22 ± -R33, and finally the 3 terms in R 00 . We have therefore exhausted 
the information available from the vacuum field equations. 

Assuming that we have a vacuum we can switch attention to the Weyl tensor, and we compute 
first 

^4 = -^0202 = [i( a 22/ 

Using [(38)1 we ma y rewrite this as 



-j - a 33 ,uu esc 2 9) + ~ia 2 3,uu esc 9] /r + 2 
* 4 = KJr + 2 , 



(46) 



(47) 



where 

A/" = \ (a 22 + ia 23 esc 9) ;U (48) 
is the Bondi news function ( Bondi et al. . 1962h . 

This is a hig hly satisfactory re sult w hich, in spite of o ur rather ad hoc chart and tetrad, mimics 
the treatment of Bondi et al.1 (|1962l) and lNewman & Until (|1962r) . Further we see that it is linear in 
the a mn and so should appear in linearized theory. Also it does not involve the Lorentz parameters 
a n and (3 n and so is tetrad-invariant (for tetrads which are asymptotically Minkowskian). The 
remainder term in in |(47)| contains some 2 terms and junk 3 terms. 

Next consider 



* 3 = R 



0120 



y { 3 2) /r 2 + 3 



where 



(2) 



2 1/2 (AT - iAf^ esc 9 + Af cot 9) + (a 4 + ia 5 )N u . 



(49) 



(50) 



Note first that the r-dependence is precisely what one would have expected from the peeling prop- 
erty. The first term in the coefficient fyf^ is linear and would have been predicted within linearized 
theory. However the second term is nonlinear for it depends on the a n which determine the in- 
finitesimal Lo rentz transformation of the NP tetrads [(30)] and |(3"Tjj This is to be expected. The NP 
tetrad used by iNewman & Until (| 19621) was chosen very specifically, while here we are considering 
a class of tetrads infinitesimally close to the Minkowski one. If we were to restrict attention to the 
subclass of tetrads where a 4 = a 5 = then our result would be consistent with linearized theory. 
On the other hand another choice of a 4 + ia 5 would give fyf^ = 0. The remainder term in |(49)| is 
junk. 

Next we find that 



^ 2 = R 



1320 



*f /r 3 + O 



where the remainder term is junk. We will return to the leading term shortly. 
We find next that 

Vi = Roii3 = y { i ) /r 4 + 5 , 



(51) 



(52) 



The coefficient ty^ contains nonlinear terms, but we are unable to determine it precisely because 
it also contains junk terms. The peeling property is still holding though. 

The peeling property would demand that \l/ — -R1313 should be 5 . However we find 



^ 4) /r 4 + O 



5j 



(53) 
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where 



^q 4) — | [( fl 22 + 033 esc 2 9) (a 22 - 2za 2 3 esc 9 - a 33 esc 2 9) + 4(b 22 - b 33 esc 2 9 - 2ib 23 csc 9)}. (54) 



But w e have not made the restrictions that were imposed by lBondi et al.l(|1962|) or lNewman & Unti 
dl962h to ensure peeling. The latter restr iction was to deman d = 0. The former restriction was 
that the "outgoing radiation condition", (|Bondi et all 1 19621) . held. In our notation this condition is 
, see |(28)1 

^22 = -033 csc 2 9, b 22 = b 33 csc 2 9, b 23 = 0. (55) 

Examining |(54)| we se e that imposing the ou tgoing radiation condition |(55)| ensures * { o ] = 0, a 
result first obtained by lValiente-Kroonl (|1999|) . With one or other condition we have 



^/r b + 6 



(56) 



(5) 

but the coefficient \I/q contains junk terms and so we cannot evaluate it. (It also contains nonlinear 
terms not predicted by linearized theory.) 

To summarize: if we impose the outgoing radiation condition |(55)| then we obtain the peeling 
property, and we can obtain explicitly the leading terms in \1>4, ^ 3 , \1> 2 , but not those for ^ and \& 
because they contain junk terms. 

We now return to the discussion of ^2 given by |(51)| The leading term coefficient is 



(3) 



:|a 00 - \ia 23 csc 3 9 + \b 2 2, u + \ (^22 + ^23 csc 9) (a 22 , u - ia 23 , u csc 9) + 

i?(a 2 3,e - 2a 2 2 l9 i) cot 9 csc 9 - \%a 22M csc 9 + \i{a 23fie - csc 2 9a 23tH ) csc 9. 



Here we have fixed the Lorentz parameters, as described earlier, and are imposing the outgoing 
radiation condition |(28)| 

Now the Bo ndi mass M B (u) can be defined by, (|Bondi et all Il962[ ). (|Newman & UntiL Il962l) . 
dStewartl[l989h . 



4irM B (u) = - lim / r 3 (^/ 2 + aX) sin9 d9 , 

r ^°° JS(u,r) 



(58) 



where the integral is over the 2-surface S(u, r) given by u = const, and r = const. Here a and A 
are NP spin coefficients given by 



a = m a 5l a = 7313 = ^ /r 2 + 3 , 

where 

<r (2) = -\[a 22 - ia 23 csc 9), 
Taking the limit in |(58)| we have 



A = n a 5m n 



7022 



A (1) /r + 



2- 



(a 22 . u + ia 2 3,n csc 1 



4nM B (u) = - [ (¥ 2 3) +a (2) A (1) ) sinfld^ 

JS(u,l) 



(59) 
(60) 

(61) 



Of course the formula [(58)] is only v alid in a spec ially chosen Bondi frame. The generalization 
to an arbitrary NP frame is discussed in Stewart (|l989|) . In the large r limit our frame differs from 
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the Bondi one by a Lorentz transformation which is close to the identity. A 2-parameter subgroup 
of the Lorentz group consists of "boosts" and "spins" 



I — *■ a L 



n 



a 2 n, 



m 



e^m, 



(62) 



where a and ip are real. Using the formulae in appendix B of [Stewart] (|1993|) it is easy to verify that 
the integrand of |(61)| is invariant under boosts and spins. Next consider a 2-parameter subgroup of 
"null rotations about I" given by 



m 



m + cl, 



n 



n + cm + cm + ccL 



where c is complex. Under such a transformation 



# 2 -»• ^ 2 + 2c^i + c 2 ^ , a -> <r + ck, 

A — > A + C7r + 2ca + c 2 (p + 2e) + c 3 k + cl a c, a + rrf'c i a . 



(63) 



(64) 



We expect c = 0\ and the NP scalars a, 7r, p and e are all Oi. Thus the integrand of |(61)| is not 
changed. We should also consider null rotations about n given by 



n 



m + cn, I — > I + cm + cm + ccn, 



so that 



^2 -> ^2 + 2C^ 3 + C 2 ^ 4 , A -> A + CI/, 

cr -> o- + cr + 2c/5 + c 2 (/i + 27) + c 3 z/ + cra a c + m a c a . 



(65) 



(66) 



Now we have taken great care to ensure that I is almost geodesic ( k — 4 ) and almost affinely 
parametrised (e + e = 3 ) and so we should only consider the transformation |(65)] where c = 3 . 
Under this restriction the integrand of |(6 1 )| is not changed. Thus the formula |(6 1 )| evaluated in our 
frame does indeed give the Bondi mass to leading order. Next note that 



(67) 



Im(&2 } + o" (2) A (1) ) = - ±a 23 esc 3 9 + |a 2 3,e esc 9 cot 9 + \a 23 ,ee esc 9- 

\a 2 2,<p esc 9 cot 9 - \a 2 2fi4> esc 9 - |a 23) ^ esc 3 6*. 

When we integrate this over the unit sphere the terms in the second line give zero since their 
contribution to the integrand is 27r-periodic in <\>. Those in the first line contribute 

\k [csc#(sin 9a2z),e\o ■ 

Now a 2 3 must scale like sin 2 9 at the end points or else the integrand is singular. It follows that the 
Bondi mass must be real, and 



AnM B (u) 



(a o + &22,« + a220 2 2,« + 023023,^ esc 2 9) sin 9 d9 < 



(68) 



Finally there is a standard result, (IBondi et all 1 19621) . (INewman & UntiLll962f) . (IStewarill989l) . 
for the rate of decrease of the Bondi mass 



a dMB ( \ 
4/T — - — (u) 

du 



\M\ 2 sin9 d9 d<p 



s(u,i) 



s(u,i) 



((fl22,n) 2 + (Q-23,m) 2 CSC 2 9) SU1 9 d9 i 



(69) 
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demonstrating the well-known result dM B /du ^ 0, the Bondi mass decreases as energy is radiated 
away, a result not deducible in linearized theory. Although |(69) | was originally derived in a special 
Bondi frame, it too holds in our approximate Bondi one, at least to leading order. We should 
empasize that although the outgoing radiation condition was used in the derivation of |(68)[ the 
mass loss formula [(69)1 holds without the need for this restriction. 
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6 Implications for numerical relativity 



At first glance the formalism we set up to carry out this study may seem to be cumbersome, but 
it has the advantage the the results can be translated back into the X a = (T, R, 9, $) chart, after 
which the theoretical chart x a = (u, r, 9, <p) can be discarded. 

We chose the x a chart so that the metric coefficients an and b\\ vanished as well as a\ 2 and a\ 3 . 
Now 



d \ fdT\ fd \ fdR\ fd 



du J _ \du J r \dT J R V du J „ \ dR 
[1 - qir- 1 ) (?-) +() 



dT 



(70) 



2, 
R 



using |(23)| In principle the function g x is arbitrary. But the vacuum field equations implied a 01 = 0, 
and then equations |(A3)| imply 

qi = | (/ioo - /in)- ( 71 ) 

The vacuum field equations imply [(38)] 

d22,u + fl33,« CSC 2 9. 

Using |(A3~)1 and [(43)] we have, to leading order 

[h 22 + 2/i 02 ,e + 2/n 2 ,e],T + [/133 + (2 V* + 2/l i3,$) csc2 9 ],t csc 2 9 = 0. (72) 

Other vacuum conditions can be handled in a similar way. 

In order to discuss the Bondi news function and Bondi mass, it is convenient to introduce some 
auxiliary functions in the numerical chart, 

W =/io3 + 

A =h 22 + 2(/i 2,e + h l2t e), 

B =h 23 + h 02 ,* + fin,* + W, e - 2W cot 9, 

C =k 22 + 2Wh 23 cot 9 csc 2 9 - 4W 2 cot 9 csc 2 9 + (k 02 + k 12 ) :B - 

(|/in + M Ve - ((|/ioo + h i)h 02 ) >e + (4Wcot 9 - /i 23 )W, e csc 2 9- 
lh 02 h lhe - h 22 h 12)Q + (h 02 + hi 2 )h 22 ,e + W/i23,e csc 2 9 - (W i0 ) 2 csc 2 9 

which should be readily available according to the assumptions in section [2l 
Then the leading term in the Bondi news function given by [(48)] becomes 

N = A iT + iB, T csc B, (74) 

whose calculation might require some sophistication, although there is no reference to the interme- 
diary x a = (u, r, 9, 0) chart. Because the news function is linear in the h a b and their derivatives, it 
could have been calculated within linearized theory. 
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In particular the formula [(68)] for the Bondi mass M B {u) translates into 



M B (T~R)=-—f [ (h ll + AA T + BE T esc 2 Q + Ct) sinQdQd^. (75) 
87T J^ =0 J e=0 V ' / 



In a similar way the formula [(69)] for the rate of change of M B at fixed large R is 



M B (T - R) 



16tt 



$=0 J 0=0 



(A iT )' z + (£ t) 2 esc 2 6 ) sin 6 dQd<5>. 



(76) 



Working in linearized theory numerical relativists often take the leading terms in A and B to 
represent the "gravitational waveforms" h + and h x . Equations |(73)| suggest that other terms are 
present, even in linearized theory. 

Again we emphasize that the intermediate x a chart does not intrude — the formulae |(74)[ |(75)| 
and |(76)| apply in the numerical X a = (T, R, 6, $) chart. Only the first can be deduced from 
linearized theory. 



Why are these formulae s o com plicated, when compared with the original papers, iBondi et al. 



( 1962 ) and Newman & Until ( 1962 ). or even the formulae in section [5]? Well the coordinates and 
tetrads of the originals were very carefully chosen to simplify the problem, and much of this paper 
has been spent building the relationship between the numerical relativist's X a = (T, R, 6, $) chart 
and the x a = (u, r, 9, <p) chart and adapted tetrad used in this paper, in which the formulae look 
almost as simple as in the original approaches. One way to a void the compl exity is to design a 
numerical approach based on Penrose's geometrical approach dPenrosd . Il963|) . but that brings in 
different problems and complexities. 
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A Computational Details 



The h and k occurring in (1 1) are given by 



h 00 = -h 00 , h m = hoi, h 02 = h 02 , h 03 = esc 2 9h 03 , h u = -h n , 
h 13 = - esc 2 Ohxa, h 22 = -h 22 , h 23 = - esc 2 9h 23 , h 33 = - esc 



and 



k 00 = -k 00 + h 2 m - h 2 m - h 2 02 - esc 2 9h 2 03 , 

k 01 = feoi - h 00 h 01 + h 01 h u + h 02 h 12 + esc 2 9h 03 h 23 , 

k 02 = k 02 - h oh 02 + h 01 h 12 + h 02 h 22 + esc 2 9h 03 h 23 , 

k 03 = esc 2 9(k 03 - h 00 h 03 + h 01 h 13 + h 02 h 23 + esc 2 9h 03 h 33 ), 

k 11 = —kn + h 01 — h\ x — h\ 2 — esc 2 9h\ 3 , 

k 12 = -k l2 + h 01 h 02 - h u h l2 - h 12 h 22 - esc 2 9h 13 h 23 , 



k 13 = esc 2 9(-k 13 + h ih 03 - h n h 13 - h X2 h 23 - esc 2 9h 13 h 3 
k 22 = —k 22 + h\ 2 — h\ 2 — h\ 2 — esc 2 9h\ 3 , 
k 23 = esc 2 9(-k 23 + h 02 h 03 - h 12 h 13 - h 22 h 23 - esc 2 9h 23 h 33 ), 

4)- 



A; 33 = esc 4 9(-k 33 + /iq 3 - h\ 3 - h 2 23 - esc 2 9h\ 



The a ab occurring in |(26)| (after imposing the conditions |(27)[ ) are given by 

000 — ^00 — 2gi, 

001 = ^00 + h i + 2m 1 - q u 
a-02 = h 02 — y 3u , 

a 3 = ^03 - z 3 ^ u sm 2 9, 

an = hoo + 2/ioi + hn + Ami, 

a>i2 = h 02 + h 12 + y 2 , 

013 = ^03 + ^13 + z 2 sin 2 9, 

a 22 = h 22 - 2y 2 e , 

a 23 = h 23 - y 2 4> - z 2fi sin 2 9, 

a 33 = h 33 — 2z 2 s. 
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(A4) 



The relation between the tetrad components c mn and the c mn is (here s = 2 1 / 2 ) 

c 00 — \cq\— Coo, 

c 01 = \cq — |c i — c 10 + |cn, 

c 02 = s(c 20 + c 30 ) - \s(c 21 + C31), 

c 03 = is[(c 2 o - c 3 o) - |s(c 2 i - c 3 i) esc (9], 

c 11 = \cqi - c n , 
c 12 = s(c 2 i + C31), 
c 13 = is(c2i — c 3 i) csc (9, 
c 20 = s(c 02 + «c 03 csc 9), 
c 21 = s(ci2 - \c m ) - is(\c m - ci 3 ) csc (9, 
c 22 = -|(c 22 + c 32 ) - |i(c 23 + c 33 ) csc (9, 
c 23 = — |i(c 2 2 - c 32 ) csc (9 + |(c 23 - c 33 ) csc 2 (9, 
c 30 = s(c 02 - ic 03 csc 9), 

C 31 = S(ci2 - |c 02 ) + «s(|c 03 - C13) csc 9, 
C 32 = -|(C 22 + C 32 ) + |«(c 23 + C33) CSC (9, 

c 33 = — |z(c 22 - c 32 ) csc 9 - |(c 23 - c 33 ) csc 2 9. 

One possible relation between the tetrad coefficients c mn and the metric coefficients a mn and 
Lorentz parameters a m is 



c 00 




«1, 


Coi 






C02 




Oi2 — 2SGC4, 


0)3 




013 + 2so;5 sin (9, 


ClO 




|aoo — |o!i, 


Cll 




a oi ~~ ~~ a i, 


C12 




ao2 — |ai2 — s(2« 2 — a 4 ), 


Cl3 




ao3 - |oi 3 + s(2a 3 - a 5 ) sin 


C20 




«2 + ^«3, 


C 2 1 




0:4 + i«5, 


C22 




^sa 2 2 — is(a 2 3 — a 6 ) csc (9, 


C23 




sctg — |isa 33 csc (9, 


C30 




a 2 — ia 3 , 


C31 




0:4 — i«5, 


C32 




\sa 22 + is(a 23 - a 6 ) csc 9, 


C33 




saQ + \isa33 csc 9, 



(A5) 



where s = 2 1 / 2 . Other representations are possible. 
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B Tetrad formalism 



Here we define the notation and summarize the results. The reader to whom this material is unfa- 
miliar should consult introductory material e.g.. IChandrasekhai dl983h chapter 1, section 7. 



In this section a, b, c, . . . are coordinate indices while a, /3, 7, . . . are tetrad indices. 
At each space-time point P we introduce a basis of vectors 

e a a , ae[0,3], ae[0,3]. 

Then the matrix 



(Bl) 



/ en 



e eo" 
ci° ex 1 



V : : 

is non-singular and we denote its inverse by e a a . Thus 



\ 



e a e p 



e n a e a , 



h- 



(B2) 



The e a a represent the dual basis of covectors. As usual chart indices are lowered (raised) using 

9a b (g ab ). 

An additional assumption made here is that 



— 9ab &a &/3 

is a constant symmetric matrix with inverse e a/3 . Thus 

The choice e a/3 = diag(l, —1,-1,-1) gives an orthonormal tetrad, but here we choose 



(B3) 



(B4) 



,a/3 



/ 1 \ 

10 

-1 

\ -1 / 



(B5) 



which gives a NP tetrad. [Newman & Penrose! (|1962|) . 
Then it is easy to see that 

e ne^ = e e a/3 e« a = e c 

so that tetrad indices are lowered (raised) using e a/3 (e af3 ). 
The Ricci rotation coefficients 7^, are defined via 

&fib;c 7A/i^C &G ci 



(B6) 



(B7) 



where the metric covariant derivative has been used, and since e a p is constant we must have 7a m ^ 
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The tetrad structure constants C 7 a/ g are defined via 

[e a ,ep] = C 7 Q/J e 7 , (B8) 
and clearly C 1 a p = C 7 r aj g]. If we let |(B8)| act on a scalar function /, note that the metric connection 



is symmetric, and use |(B7)[ then it is easy to see that 
which implies 

'yXfiu ^(^^M Cxfiu C ' flux) ■ 



(B9) 



(BIO) 



It is important to realise that the C x ^ u do not involve the connection. For \(B9)\ |(B7)| and the 
fact that the connection is symmetric means that 



r x - p x Jp a p b — p b p a \ 

^ ii.v c a,oV°£i t-i/ tyx t-i/ )■ 



Next the Ricci identity applied to e a b gives 

''■(iffy 6 To/37, a&5 Ta/3<5,a^7 



and finally 



R 



cry 



a(3-yS • 



(Bll) 



(B12) 



(B13) 



Note that throughout the paper we use these forms for the curvature tensors. Thus i? i2 means R a p 
with a = 1 and (3 = 2, which is not the same as R a b with a = 1 and 6 = 2. 



Our algorithm starts from the sets {e Q a } and {e a a }. We compute C a p 1 from |(B lT)] and of 
course C 7Q/ 3 = e 7< sC£g. Next we compute 7 Q/37 from |(B 10j| and finally the curvature tensors from 



(B 12)| and |(B13)[ Although this looks ponderous it can easily be automated using a computer 
algebra system. 
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